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Abstract. The production of electron pairs with pt between 1 and 5 GeV/c 
q . and to < 300 MeV has been measured at mid-rapidity in ^/s NN = 200 GeV 

^ ' P+P and Au+Au collisions by the PHENIX experiment at RHIC. A significant 

Ch ■ excess above the hadronic background was observed in both p + p and Au+Au 

collisions. Treating the excess as internal conversion of direct photons, the 
direct photon yield in Au+Au was found to be enhanced compared to the 
^ ' binary-scaled p + p yield. The enhancement is consistent with an exponential 

inverse slope of 221 + 23+18 MeV and predictions from hydrodynamical models 
with initial temperature between 300 and 600 MeV at formation times of 0.6- 
0.15 fm/c. 
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1. Introduction 

A multitude of results from the Relativistic Heavy Ion Collider (RHIC) at Brookhaven 
National Laboratory (BNL) indicates the creation of a high-density, thermalized 
medium in ultra-relativistic collisions of heavy ions [Q]. Such a medium is expected 
to radiate thermal photons which, once produced, leave the medium unscathed. 
Thermal photons from the partonic phase of the collision are predicted to dominate 
the direct photon spectrum in the transverse momentum (jpr) range of 1-3 GeV/c 
as illustrated in Fig. |T] [ 2J . However, here direct photons are submerged in a back- 
ground of hadronic decay photons, mainly from the ir° and 77. This background 
constitutes a major experimental challenge in the conventional, calorimeter-based 
measurement. It can be overcome, though, by measuring low- mass electron pairs in 
a mass range where electron pairs from the ir° Dalitz decay do not contribute [ [3] . 
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q [GeV] mass dependence for direct photons as 

Fig. 1. Direct photon sources at RHIC we ll as for tt° and r\ Dalitz pairs. 
[®. 



The electron-pair yield above the hadronic background can be treated as internal 
conversion of direct photons [ [4] . 

2. Internal Conversion Method 

Any source of high-energy photons also emits virtual photons with very low mass 
[2]. Those virtual photons then convert to low-mass e + -e~ pairs, which can be 
measured (Internal Conversion Method). The pair yield per direct photon falls with 
the pair mass as: 

am Sir m V m z V / 

Here a is the fine structure constant, m e and m are the masses of the electron and 
the e + -e _ pair, respectively, and S is a process dependent factor that goes to 1 as 
to — > or m <C pt- Equation [1] also describes the relation between photons from 
hadron decays (e.g. 7r°, r\ — > 77, and uj — > jit ) and e + -e~ pairs from Dalitz decays 
(it , rj — > e + -e~7 and u> — > e + -e~ir a ). For 7r° and 77, the factor S is given by S = 
|F(m 2 )| 2 (l — -jfiY [IS], where Mh is the meson mass and F(m 2 ) is the form factor. 
Figure [5] illustrates the pair yield mass dependence for direct photons as well as for 
w° and r\ Dalitz pairs. The cut-off of the 7r° pairs at the 7r° mass can be exploited 
to increase the signal-to-background ratio from 10%, where it is comparable to the 
systematic uncertainty and therefore not significant, to 50%, making a significant 
measurement possible. Since the measurement at low px is systematics limited, the 
simultaneous reduction in statistical significance is an acceptable trade-o 



lr There are about 0.001 virtual photons with m ee > M n o for every real photon. 
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3. Data Set And Backgrounds 

The analysis is based on two data sets: Au+Au at y/s NN — 200 GeV acquired 
in 2004 consisting of 0.8 billion minimum bias events (4.9 pb _1 p + p equivalent); 
p + p at the same cms energy acquired in 2005 with 2.25 pb _1 . Charged tracks 
were measured with the Drift Chamber and Pad Chamber of the PHENIX [ [6] 
Central Arms covering \rj\ < 0.35 and A<fr = 2 x ir/2 and identified as electrons 
with the Ring Imaging Ccrcnkov Detector (RICH) and Electromagnetic Calorimeter 
(EMCal). Material conversion pairs were removed by a cut on the orientation of 
the pair plane with respect to the magnetic field. Combinatorial background was 
removed by an event mixing technique, accurate to 0.25% systematic uncertainty 
in Au+Au. 

There is additional correlated background from two sources: cross pairs with 
one electron/positron from either virtual photon in a double-Dalitz decay; jet pairs 
from two different Dalitz decays within the same jet or from back-to-back jets. 
These contributions can be well understood in a Monte Carlo calculation and have 
been subtracted. 

4. Signal Extraction 

After subtraction of the combinatorial background and the cross and jet pairs, the 
pair mass spectrum is compared to a "cocktail" of known hadronic sources [ [7J H] ■ 




m,.,- (GeV/c J ) nv, (GeV/c 2 ) 

Fig. 3. Pair mass spectra for data (points) and hadronic cocktail (lines) in the 
PHENIX acceptance for different px intervals for p + p (left) and Au+Au (right) [ 
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Figure [3] shows this comparison for different pt intervals for both p+p and Au+Au. 
The cocktail is normalized to the data for m < 30 MeV, where the ir° Dalitz decay 
dominates the yield. The "knee" at 100 MeV comes from the 7r° cut-off leading to 
the 80% background reduction mentioned above. 

In p + p, the pair yield is consistent with the hadronic background for the 
lowest pt interval. At higher pt a small excess is visible for m > m T « . In Au+Au a 
much larger excess appears at all pt, indicating an enhanced production of virtual 
photonfo 

To quantify the direct photon fraction, the mass spectrum is fit with a two- 
component function f(m ee ) = (1 — r)f c (m ee ) + rfdir(mee) as illustrated in Fig. [4] 
Here f c (m ee ) is the shape of the cocktail mass distribution (shown in Fig. [3]), 
fdirijn-ee) is the expected shape of the direct photon internal conversion, and r 
is the fit parameter. Both / c (m ee ) and fdi r {rn ee ) are separately normalized to the 
data for m ee < 30 MeV/c 2 , where their shapes are nearly identical. This preserves 
the meaning of r as the real direct photon fraction. From the agreement between 




data and fit it can be concluded that the data matches the expected shape for direct 
photons. 

The so extracted direct photon fraction, r, is plotted as a function of px in 
Fig. [5l While in case of p + p the direct photon fraction is consistent with pQCD, 
for Au+Au r is enhanced above pQCD. 

As this measurement is based on shape differences to extract the direct photon 
fraction, the rj/w ratio is the largest source of systematic uncertainty. This results 



2 This excess is in a different kinematic region (higher pt, lower mass) than the low-mass 
enhancement reported in [[7], which is expected to be dominated by the hadronic gas phase. 
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P t (GeV/c) p 7 (GeV/c) 

Fig. 5. Direct photon fraction, r, as a function of pr for p + p and Au+Au events 
compared to a pQCD calculation for three different scales. In the case of Au+Au 
the pQCD result is scaled by the nuclear overlap function, Taa [IS]- 



in a 7% uncertainty in p + p and 17% in Au+Au. Other sources contribute only a 
few percent as the cocktail is normalized to the data and no absolute normalization 
is required. 

In the next step, r is converted into the direct photon yield as dN dlr (pT) = 
r x dN mcl (pt ) . The inclusive photon yield for each px bin is determined by dN" Lcl = 
N data x (djv^/N e c e ), where N d f a and N° e are the measured and the absolutely 
normalized cocktail e + -e~ pair yields, respectively, both for m ee < 30 MeV/c 2 ; 
and dN^ is the yield of photons from the cocktail. Here we use the fact that the 
ratio of the photon yield to the e + -e~ pair yield for m ee < 30 MeV/c 2 calculated 
from Eq. [1] is the same within a few percent for any photon source. The systematic 
uncertainty of ^ mcl is 14% from the e + -e~ pair acceptance. 

5. Results 

Figure [5] shows the invariant yield of direct photons as a function of px for p+p and 
three different centrality classes in Au+Au. The p + p data is again compared to 
the pQCD calculation. The calculation is consistent with the data for px > 2 GeV. 
A modified power law, A pp (l +p 2 -/^) _ ™i f^s the data over the entire pr range. 
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5.1. Significance of the Modified Power Law 

The modified power law appears to yield an at least as good description of the data 
at low pt as the pQCD calculation. What is the significance of this observation? 
It is obvious that the power-law behavior of hard scattering has to break down for 
Pt — * 0. For hadrons, soft production sets in with an exponential slope. This is 
illustrated in Fig. [51 which shows a parameterization of 7r° production in p + p [ 




P T [GeV/c] P T (GeV/c) 

Fig. 6. Parameterization of -k produc- Fig. 7. Fit of direct photon cross sec- 
tion in p + p [ [9] and its various contri- tion with both a power law and a mod- 
butions (for details see text). ified power law. 



9 . The parameterization is the sum of a power law and an exponential with a 
Woods-Saxon transition between the two. The two functions are shown separately. 
The exponential dies out at high px while the power law diverges at low p?. The 
hard-scattering contribution can be understood as the difference between the data 
parameterization and the exponential contribution. One can see that it flattens out 
towards low px, deviating from the power law. This flattening corresponds to an 
onset of hard scattering as illustrated by the ratio of the hard scattering contribution 
to the power law. 

This onset is not directly observable for hadrons as the production at low pr 
is dominated by soft physics. Direct photons, however, are only produced in hard 
scatterings. This makes the onset of hard scattering at low px directly measurable. 
The onset has also been measured in Drell-Yan production [110). 

To evaluate the statistical significance of the onset, the PHENIX data were 
fitted with both a power law and a modified power law. As shown in Fig. [JJ the 
modified power law yields a smaller reduced x 2 than the pure power law. 

5.2. Au+Au Enhancement 

For 1 < pt < 3 GeV, the direct photon invariant yield in Au+Au is enhanced above 
the T/iA-scaled modified power law that was fit to the p + p data (Fig. [5]). It can be 
well described, however, if an exponential is added. The resulting fit yields negative 
inverse exponential slopes of about 220 Me\ff|. If the medium were static, T could 

3 If the p + p data are fit with a pure power law, T increases by 24 MeV in central events. 
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Fig. 8. Invariant yield of direct photons as a function of pr for p + p and three 
different centrality classes in Au+Au (solid symbols). The result of an earlier EMCal 
measurement is also shown (open symbols) [ lilt 112) . The p + p data is compared 
to a pQCD calculation shown as three lines for different scales. The dashed line is 
a fit of a modified power law to the p + p data. The Au+Au data are compared to 
the T^-scalcd p +p fit (dashed line). The solid line shows the result of a fit where 
an exponential is added to the T^-scaled p + p fit [ [3] . 

be interpreted as its temperature. For a more realistic temperature estimate, the 
data is compared to hydrodynamical models. Models that fit the data assume initial 
temperatures, Tj, between 300 and 600 MeV at formation times, To, between 0.6 
and 0.15 fm/c, where the temperature and the formation time are anti-correlated. 
The 221 ± 23(stat.) ± 18(sys.)MeV obtained from the data alone serves as a lower 
limit. Even the lower limit, though, lies above the critical temperature predicted 
by lattice calculations. 

6. Conclusion 

An excess of low-mass e + —e~ pairs above the hadronic background was observed at 
intermediate px (1 GeV < pr < 5 GeV) for both p + p and Au+Au collisions. The 
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excess can be understood as the internal conversion of direct photons and shows the 
expected 1/m dependence. In p + p. the direct photon yield is consistent with the 
result of a pQCD calculation. In Au+Au, the yield is much larger. It is enhanced 
above the binary-scaled p + p yield, represented by a modified power law fit to 
the p + p data and scaled by Taa- The Au+Au yield can be well described if an 
exponential is added to the binary-scaled p + p fit. The negative inverse slope of 
the exponential is 221 ± 23 (stat.) ± 18 (sys.) MeV in central Au+Au. For a 
static medium, this could be interpreted as the temperature. For an expanding 
medium, it serves as a lower limit of the temperature. It is well above the critical 
temperature of about 170 MeV. Hydrodynamical models that fit the data assume 
initial temperatures, Ti, between 300 and 600 MeV at formation times, r , between 
0.15 and 0.6 fm/c. Together with the earlier WA98 measurement [ 113] . this result can 
be interpreted as the first experimental evidence that strongly interacting matter 
can exceed the Hagedorn temperature of 170 MeV [[Ti], 
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